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Abstract 

The (2+l)-dimension Klein-Gordon generalised equation is numerically solved through the fi- 
nite differences method. Only the sine-Gordon case is focused: kink and antikink solutions are 
obtained in cartesian coordinates and evidence of interaction in kink-kink collision is looked for 
in propagation velocity. Then the change of shape in light bullet solutions is quantified during 
propagation and in head-on collision. Lastly, the robustness of light bullets is verified in head-on 
collisions with kink, antikink, standing kink and standing breather. A 30°-collision between a 
light bullet and a standing kink is simulated as well. 

Key words: sine-Gordon equation, soliton, light bullet, kink, standing breatlier 
PACS: 02.30.Jr, 02.70.Bf, 03.75.Lm, 42.81. Dp 



1. Introduction 

The Klein-Gordon generalised equation is written as 

where F is an arbitrary function of (j). The well-known wave equation is obtained with 
= 0, while F{4>) = (^) </> leads to Klein-Gordon equation that describes a particle 
of mass m and spin in relativistic quantum mechanics. Both cases are linear and only 
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the latter is dispersive [1]. An important nonlinear case is characterised by F{(f)) = 



2 I 1 9^(6 fm,c\'^ 



19 Equation (2) is the so-called sine-Gordon equation and allows soliton-like solutions. 

20 This type of solutions is increasingly important in the description of (at least partly) 

21 particle-like objects [2,3], which is easy to understand since solitons present very localised 

22 momentum and their shape consistency results in effective transportation of energy. For 

23 instance, the 7r-mesons are quark-antiquark bound states and some works (e.g. [3]) have 

24 been able to estimate pion mass values in good agreement with the experiment by de- 

25 scribing pions as breather-like solutions of sine-Gordon equation. In other words, the 

26 quark may be seen as a soliton (e.g. kink) and the antiquark as an antisoliton (e.g an- 

27 tikink) so that tt is a bound state of the two. Another example is the Josephson junction, 

28 where the phase difference between the electronic densities of the superconductors may 

29 be modeled by sine-Gordon equation which previews the vortex dynamics in type II su- 

30 perconductors [4,5]. The sine- Gordon equation finds also application in classical systems 

31 such as coupled pendula [6] and, finally, light bullet solutions may be identified with 

32 optical pulses propagating in different media [7,8]. 

33 In this paper, section 2 briefly describes the numerical method used, while section 3 

34 presents the results of the simulations and is organised as follows. In 3.1 the propagation 

35 and collision of kink-likc solitons is analysed and in 3.2 an analogous study is performed 

36 for light bullet solutions. Section 3.3 characterises the simulation of the collision between 

37 light bullets and other sine-Gordon solutions and, finally, in section 4 the most important 

38 remarks are drawn. 



39 2. Numerical method 

In bidimcnsional cartesian coordinates, the (2+l)-dimension Klein-Gordon generalised 
equation becomes -I- — ^^^r = F{<P)- With the normalisation substitutions x = 
—x', y — —y' and t = -^t' and loosing the primes, the following equation results: 

^'* + S-S = (^)^W (3) 



dx'^ dy"^ dt"^ 



mc 



40 In order to numerically solve (3), the finite differences method is implemented in Math- 

41 ematica 5.0 @. A rectangular domain {{t, x, y) : te [ti, tnj A xe [xi,Xn^] A ye [yi,yny] } is 

42 chosen and an UfUx- n^-point space-time grid is created with uniform steps At = *"* 



43 Ax = ^"^Z^^ and Ay — ^f-zj'- Second-order derivatives are approximately computed 

44 through centered second-order differences. 

45 Initial conditions are given by the definition of (j) and ^ at time ti . As for boundary 

46 conditions, for all simulations presented here fixed conditions are used along y-direction 

47 and periodic ones along x-direction. 

48 One must take into account the convergence conditions of the method, '^^l < | and 

49 — 5' ^^^^ in the case of space-periodic solutions, the spatial grid's resolution has to 

50 be adequate: a solution with spatial periodicity through (a,b)-direction and with period 
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51 
52 
53 
54 
55 
56 



S should be represented with some - say, 3 - points per period, that is, T = > 3. 

Otherwise, meaningless results may be obtained. 

The whole numerical approach reffers to Klein-Gordon generalised equation. However, 
in this paper, only the sine-Gordon case is studied; thus, from now on F{(f>) = (^) sine/). 
Future developments may consider any other case using the same numerical method with 
the corresponding definition of F. 



3. Results 



3.1. Kink-like solitons 



Theoretical approaches to (l+l)-dimension sine-Gordon equation [9,10] allow one to 
write an 1-soliton analytical solution in the (2+l)-dimension case as 



(j)(t, ~ Aarctg j Sign{m) 



(4) 
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where (3 ■ 



/ 2 1 — — * 

^™ ~ , |m| > 1, (i is the normalised propagation direction and tq the position 
of the center point at i = 0. 

The case m > 1 corresponds to the so-called fcmfc, an ascending step propagating 
along d with velocity /?. An example is presented on figure 1. Analogously, m < — 1 
corresponds to an antikink, a descending step propagating along d with velocity —f3 (see 
figure 2). Both kink and antikink behave like solitons since they maintain their shape 
while propagating. Particular solutions are the cases m = 1 {standing kink) and m = — 1 
{standing antikink) which are stationary. 




Fig. 1. Propagation of a kink with m = 2, d=ex and ro=0 within an 800 ■ 401 ■ 2-point grid defined in 
{(t, X, y) : te [0, 100] A xe [—40, 40] A ye [—40, 40]}. In the sequence of shots time flows from left to right. 



67 3.1.1. Propagation 

6B The dependence on the m parameter of the kink solution is studied. Regardless of 

69 the value of to, a nonlinear trail dued to numerical reasons is always present in the 

70 propagation of the kink. As m rises, the trail becomes more significant, which may be 

71 explained by the steepening of the step. Indeed, (j> oc arctg V'^-fi'^ ^ — arctg (e'™'^) 

72 I d=ex,ro—0 ). For m = 7 or higher, the kink is destabilised and no proper propagation 
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Fig. 2. Propagation of an antikink with m = —2, d=ex and ro=0 within an 800 ■ 401 ■ 2-point grid 
defined in {(t,x,y) : te [0, 100] A [-40, 40] Aye [-40,40]}. 



is achieved. Another important m-dependent feature is the propagation velocity. The 

theoretical motion is given by x = xo + (3t (if d=ex^ so, the velocity should be constant 

and equal to f3. To determine this velocity numerically, several pairs of time and position 
of the center point are registered and a linear fit is applied using Origin 5. @- an example 

is presented in figure 3. The results for m = {f .5, 2.0, 2.5, 3.0} {using d=e-x-,'rQ=Q^ are 
presented in table 1. It is obvious that greater values of m lead to less precision and 
exactness on the determination of the velocity. Taking into account the last column of 
table 1, the most favourable situation occurs when m = 2.0. Therefore, the next sections 
will use this kind of kinks (and the corresponding antikinks with m = —2.0) as they are 
less disturbed by numerical errors. IVLoreover, the fitted values of xq are close to 0, which 
means the law x — xi^ + /3ct is being followed {^0= ^ ■ 

It is interesting to note that all velocities in table 1 are below 1 (c = 1), which is the 
typical velocity of a wave in the case = 0. As m ^ oo, the velocity of the kink 

tends to 1 since (3^1. 




Fig. 3. Linear fit of pairs {time, position} of the center point of a kink with m = 2, d=£x and ro=0 
propagating within an 800 ■ 401 ■ 2-point grid defined in {(t, x, y) : te [0, 100] A xe [-40, 40] A ye [-40, 40]}. 



m[] 






^(teo) [] 


-5/3 [%] 


1.5 


-0.093 ± 0.000 


0.743 ± 0.000 


0.745 


0.28 


2.0 


-0.108 ± 0.000 


0.864 ± 0.000 


0.866 


0.26 


2.5 


-0.081 ± 0.028 


0.908 ± 0.001 


0.917 


0.97 


3.0 


0.017 ± 0.039 


0.922 ± 0.002 


0.943 


2.25 



Table 1 ^ ^ 

Determination of the propagation velocity of kinks with m = {1.5, 2.0, 2.5, 3.0}, d =6^ and ■ro= prop- 
agating within an 800 ■ 401 ■ 2-point grid defined in {{t, x, y) : te [0, 100] A xt [-40, 40] A ye [-40, 40]}. The 

last column shows the relative difference between ^{"i"") and /3(*'^°): 5^ = -I p(t<io) ' 

87 3.1.2. Collision 

The collision between kink- like solitons is achieved by superposing them. However, in 
general, this procedure is not legitimate since the sine-Gordon equation (2) is nonlinear. 
The condition for the coexistence of a pair of legitimate kink-like solitons is 

sin{(j)i + (j)2) — sin(j)i -f sin4>2 (5) 

8B where (pi and 02 £^re solutions of (2). Thus, two kink- like solitons may be superposed if 

89 their steps are sufficiently separated, because in that case on each grid point one of the 

90 two solutions is a multiple of 27r, which means (5) is verified. In this way, the kink-kink 

91 collision, shown in figure 4, can be easily started and studied. Only the collision of kinks 

92 with symmetric velocity and same |m| is simulated. Future works may take into account 

93 other situations. 




Fig. 4. Kink-kink collision. Both kinks have m = 2 and d=ex\ the kink propagating left-to-right (K— >) 
has ro= {— 10, 0), while the one propagating right-to-left (K<— ) has ro= (10,0). An 800 ■ 401 ■ 2-point 
grid is defined in {{t,x,y) : te [0, 100] A xe [-40,40] A ye [-40,40]}. 

94 Figure 4 shows that the kinks emerge intact (shape-wise) from collisions, which is 

95 typical of a soliton's behaviour. However, as they are nonlinear objects, they interact with 

96 each other during collisions. Evidence of this interaction is looked for in solitons velocity. 

97 The determination of the velocity follows the procedure explained in section 3.1.1, but 

98 picking pairs of time and position only after the interaction. The results obtained for 

99 each of the kinks are presented in table 2. Comparing these results with those referring 

100 to propagation only (table 1, line with m = 2.0), one understands that the value of /3 does 

101 not change significantly, while xq is now clearly above 0. This means that the interaction 

102 speeds up each one of the kinks, but does not alter their velocity afterwards. 
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0.703 ± 0.000 


0.864 ± 0.000 




0.703 ± 0.000 


0.864 ± 0.000 



Table 2 

Determination of the propagation velocity of the kinks in the kink-kink collision referenced by figure 4. 

103 3.2. Light bullets 

104 The light bullet [7,8] represents a well-localised two-dimensional moving pulse. Its pre- 

105 viewed space-time dependence is 



4> (t, 7^ = sin (^"/ d ■ (j^ - n,^ + Lot^ ■ 

-3^([d-(-"^o)] +p-(r-ro)]') 



(6) 



106 where uj = ^1 + 7^ and u is a normalised direction perpendicular to d ■ 



107 3.2.1. Propagation 

The evolution of a single light bullet is represented in figure 5. There is no significant 
modification of shape - light bullets behave indeed as solitons. Nevertheless, a nonlinear 
trail is formed and the cr-value of the envelope function changes during propagation. To 
quantify the former, one compares the integral of (j) in the trail region to that in the pulse 
region: 

[ , (hdxdy 

The higher e, the more important the trail is in comparison with the soliton. It is im- 
portant to note that the trail is dued not only to numerical reasons but also to the fact 
that (6) is not an analytical solution of (2), unlike the case of the kink-like solitons. The 
latter effect is evaluated fitting the numerical values of (j) after propagation to 



^(p-(r-al)]'+[n.(7-a3)]') 



oo + Oie 

108 where sometimes not all parameters are free to vary for convergence purposes. 




800 • 151 • 15-point grid is defined in {(t, x, y) : te [0, 100] A xe [-30, 30] A ye [-30, 30]}. The images are the 
projection of the solution in the xz plane. In this case, F = 7.85. 
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109 One can now study e and a2 as functions of 7 and a and then choose the most favourable 

110 pair {7, cr} to perform collisions between light bullets. Tables 3 and 4 present the results 

111 of these simulations. Although the case {7, a} = {1.0, 5.0} is the one with less significant 

112 trail, it leads to a strong deviation of a after propagation. The pair of parameters {7, a} = 

113 {2.0, 2.5} seems to be the one in which less numerical errors occur; therefore, these are 

114 the values used in the next sections. 





1.0 


2.0 


3.0 


1.0 


0.0402 


0.0439 


0.0447 


2.5 


0.0448 


0.0430 


0.0442 


5.0 


0.0380 


0.0487 


0.0461 


r[] 


15.71 


7.85 


5.23 



Table 3 

Values of e for light bullets with 7 = {1.0,2.0,3.0} and a = {1.0,2.5,5.0} propagating within an 800 • 
151 • 15-point grid defined in {{t,x,y) : te [0, 100] Ai;e[— 30, 30] A j/e [— 30, 30]}. F-values are also shown 
for each 7 used. 



<^ D \ 7 [] 


1.0 


2.0 


3.0 


1.0 


2.688 


1.156 


2.067 


2.5 


2.704 


2.468 


2.573 


5.0 


5.907 


4.596 


4.639 


r[] 


15.71 


7.85 


5.23 



Table 4 

Values of the fitting parameter 02 for light bullets with 7 = {1.0,2.0,3.0} and a = {1.0,2.5,5.0} 
propagating within an 800 •151- 15-point grid defined in {{t, x, y) : te [0, 100] A xe [—30, 30] A ye [—30, 30]}. 
F-values are also shown for each 7 used. 



115 3.2.2. Collision 

116 Only the collision between light bullets of {7, a} = {2.0, 2.5} is simulated. Other cases 

117 may be delt with in future approaches. The numerical method is started by superposing 

118 light bullets according to (5). Figure 6 shows a head-on collision. The solitons interact 

119 with each other and emerge esscintially intact and with no change in propagation direc- 

120 tion. However, the cr-value of the envelope function of each bullet does not behave as in 

121 the propagation case. In fact, the value of the fitting parameter 02 after collision is 2.09. 

122 It seems that head-on collisions lead to a shrink of the soliton, but this conclusion must 

123 be tested by further and more complete simulations. 

124 3.3. Collision between light bullets and other solutions 

125 In order to study the interaction of light bullets with kink-like solutions, kink-, antikink- 

126 and standing kink-light bullet head-on collisions are set up and shown, respectively, in 

127 figures 7, 8 and 9. In all three cases, the light bullet seems to emerge intact from the 

128 interaction and continues its motion with no visible change. The velocities of the light 
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Fig. 6. Head-on collision between two light bullets with {7,0"} = {2.0,2.5}, roi= (—15,0) and 
r^2= (15,0). An 800 ■ 151 ■ 15-point grid is defined in {(t, x, j/) : te [0, 100] A xe [-30,30] A ye [-30, 30]}. 
The images are the projection of the solution in the xz plane. 

129 bullet in the propagation case (section 3.2.1) and in kink- light bullet head-on collision 

130 arc roughly determined and no significant difference is noticed. 

131 Moreover, a 30°-collision between a light bullet and a standing kink is simulated - see 

132 figure 10. As in the previous scenarios, no modification in the light bullet propagation 

133 direction is detected, which evidences the robustness of these objects. 




Fig. 7. Kink-Ught bullet head-on collision. The light bullet has {7,0-} = {2.0,2.5}, roi= (-15,0) 

and di=ex, while the kink presents m = 2, ro2=0 and ^2= — £x- An 800 • 251 • 31-point 
grid is defined in {(t, a;, y) : [0, 100] A xe [— 30, 30] A j/e [— 35, 35]} and this domain is zoomed into 
{(t, X, y) : te [0, 100] A xe [—20, 15] A ye [—15, 15]}. In the sequence of shots time flows from left to right 
and downwards. In this case, P = 13.09. 

Another collision taken into account is the one between a light bullet and a standing 
breather, shown in figure 11. The latter object is an analytical solution of (2) which is 
oscillatory and may be interpreted as a bound state between a kink and an antikink. Its 
behaviour is previewed analytically in the (l-l-l)-dimension case [9] and can be extended 
to 

, / -►x , / TO sinujt \ 

(t, r ) = Aarctg / ^ ^xn (7) 

^ ^ coshirn d -{r -ro\\ j 

134 where = ^ = \/l — m? and |to| < 1. In this collision, the light bullet still emerges 

135 intact, although the standing breather is completely ruined after the interaction. In other 
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Fig. 8. Antikink-light bullet head-on collision. The light bullet has {7,0-} = {2.0,2.5}, roi= (15,0) 

and di= — e^, while the antikink presents m = —2, ro2=0 and d2= — e^. An 800 ■ 251 ■ 31-point 
grid is defined in {(t, x, y) : [0, 100] A xe [— 30, 30] A j/e [— 35, 35]} and this domain is zoomed into 
{it,x,y) : te [0, 100] A [-15, 20] A ye [-15, 15]}. In this case, T = 13.09. 




grid is defined in {(t, x, y) : te [0, 100] A xe [— 30, 30] A j/e [— 35, 35]} and this domain is zoomed into 
{{t,x,y) : te [0, 100] A xt [-20, 15] A ye [-15, 15]}. In this case, T = 13.09. 



136 words, the bound state kink- antikink is destroyed by the hght bullet, even though this 

137 does not destroy isolated kinks nor antikinks. 

13B Figures 7, 8, 9, 10 and 11 all together allow one to confirm the robustness of light 

139 bullets observed in 3.2. Indeed, these solutions pass through analytical solutions (4) and 

140 (7) and remain unchanged. This property may have interesting consequences, specially 

141 if one identifies light bullets with optical pulses propagating in different media [8] . 

142 4. Final remarks 

143 The analysis of kink and antikink solitons revealed an obvious shape consistency during 

144 propagation, which is typical of solitons behaviour. JVIoreover, an m-dependent study of 

145 the velocity in the kink case was carried out and it followed the theoretical solution (4), 

146 as expected. 
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Fig. 10. Standing kink-light bullet 30°-collision. The Ught bullet has {7,0-} = {2.0,2.5}, 

roi= (— 15cos30°, — 15sin30°) and di= cos30° +sm30° Sy, while the standing kink (m = 1) presents 

r^2=~0 and d2=e*x- An 800-61-61-point grid is defined in {it,x,y) : te [0, 100] A xe [-30,30] A ye [-30, 30]} 
and this domain is zoomed into {{t,x,y) : te [0, 100] A xe [-20, 20] A ye [-20, 20]}. In this case, T = 3.14. 




Fig. 11. Standing breather-light bullet head-on collision. The light bullet has {7,0"} = {2.0,2.5}, 

roi= (—10,0) and di=ex, while the standing breather presents m = 0.8, ro2=0 and d2=ex. An 
800 ■ 251 ■ 31-point grid is defined in {(t, x, y) : te [0, 100] A xe [-30, 30] A ye [-35, 35]} and this domain is 
zoomed into {{t,x,y) : te [0, 100] A xe [-25,25] A ye [-15, 15]}. In this case, V = 13.09. 

147 As for kink-kink collision, one verified that kinks do not alter their shape, propagation 

14B direction nor velocity after head-on collisions with each other. Nevertheless, the numerical 

149 value of was found to be significantly above for both kinks, which proves that 

150 nonlinear interaction has taken place and that it momentarily speeded up both solitons. 

151 In the light bullet case, the propagation was characterised as a function of the param- 

152 eters 7 and a and the most favourable case was identified: {7. ct} = {2.0, 2.5}. With this 
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153 configuration, a head-on collision was set up and nonlinear effects in the shape of the light 

154 bullets were measured (through 02). The shape consistency of the light bullets noticed in 

155 the latter collision was also seen in head-on collisions with kink, antikink, standing kink 

156 and standing breather. 

157 Moreover, a 30°-collision between a light bullet and a standing kink was set up as 
15B well and it reinforced the idea that light bullets and kinks are indeed robust objects. We 

159 also found that the direction of the propagation of the light bullet is maintained after it 

160 collides with the kink. 

161 Future developments may study the dependence on 7 and a of light bullet collisions 

162 and verify if these lead to a shrink of the light bullets or not. Another interesting point is 

163 to understand why kink, antikink and standing kink survive to the collision with a light 

164 bullet while the standing breather does not. 



165 
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